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Abstract 

Observations of blazar flaring states reveal remarkably different vari- 
ability time scales. Especially rapid flares with flux doubling time scales 
of the order of minutes have been puzzling for quite some time. Many 
modeling attempts use the well known linear relations for the cooling and 
emission processes in the jet in a steady-state scenario, albeit the obvi- 
ous strongly time-dependent nature of flares. Due to the feedback of self- 
produced radiation with additional scattering by relativistic electrons, the 
synchrotron-self Compton (SSC) effect is inherently time-dependent. Al- 
though this feedback is usually implemented in numerical treatments, only 
recently an analytical analysis of the effects of this nonlinear behaviour has 
been performed. Here, we report our results concerning the effect of the 
time-dependent SSC on the spectral energy distribution (SED) of blazars. 
We calculated analytically the synchrotron and the SSC component, giv- 
ing remarkably different spectral features compared to the standard lin- 
ear approach. Adding an external photon field to the original setting, we 
could implement quite easily the effect of an additional external Compton 
(EC) cooling, since such strong external photon fields are observed in flat 
spectrum radio quasars (FSRQ), a subclass of blazars. Calculating the 
resulting flux due to the EC cooling, we were able to show that the result- 
ing inverse Compton component strongly depends on the free parameters, 
and that SSC could potentially have a strong effect in FSRQs, contrary 
to what is usually assumed. 

1 Introduction 

Blazars, as a class of Active Galactic Nuclei (Urry & Padovanni 1995), are 
characterized by strong variability in all wave-bands. Additionally, a broad 
range of variability time scales is observed, ranging from minutes, as in the 
cases of PKS 2155-304 (Aharonian et al. 2007) and PKS 1222+216 (Tavecchio 
et al. 2011), to months (e.g. Abdo et al. 2010). Especially the very short 
variability time scales are puzzling, since they imply emission regions to be 
much smaller than the event horizon of the black hole, which should be a lower 
limit on the jet width. 

Under such circumstances it seems quite surprising that most modeling at- 
tempts still use the usual linear approach in a steady-state scenario (Ghisellini 



et al. 2009). Especially rapid flares are far from steady-state, and need a 
real time-dependent treatment. Numerical studies normally incorporate time- 
dependent features by default. There is one process, where this is particularly 
important: the synchrotron-self Compton effect (SSC, Jones, O'Dcll & Stein 
1974). Under the influence of a magnetic field relativistic electrons generate 
synchrotron radiation, which is afterwards up-scattered to very high energies 
by the same electrons. Since the electrons naturally lose energy during this 
process the energy density in the synchrotron field will decrease over time, and 
so will the efficiency of the SSC process. Since numerical treatments calculate 
the synchrotron and SSC radiation after every time step, this time-dependency 
is implemented almost unintentionally. 

Only recently Schlickeiser (2009) calculated analytically the fully time-dependent 
cooling term for the SSC radiation, which depends on an integral over the elec- 
tron distribution function itself. The differential equation (Kardashev 1962) 

^^ L -^[\itotH 1 ,t)} = s{i,t), (i) 

describing the behaviour of the electron distribution n(^,t) under the influence 
of cooling described by 7 tot and the source £(7, t), then becomes nonlinear. 

Due to the time dependency the rate of the SSC cooling decreases over time. 
This implies that, although the SSC cooling might be stronger in the beginning 
than the linear synchrotron cooling, the synchrotron cooling eventually becomes 
stronger than the SSC cooling. This can be incorporated by using both cooling 
terms in the differential equation Q , giving 



l7*ot| = Hsynl + |7W| - A)7 2 + A Ql 2 I d 77 ^(7,t), (2) 



where 7 is the Lorentz factor of the electrons, Do = 1.256 x 10 _9 6 2 s _1 , Aq — 
1.15 x 10~ 18 i?i56 2 cm 3 s -1 , with a magnetic field of strength B = frGauss, and a 
radius of the emission blob R = 10 15 i?i5cm. 

In section 2 we will present the solution to the differential equation , which 
depends on the injection parameter a, which is the ratio of SSC to synchrotron 
cooling at time of injection. We will also briefly discuss the implications of this 
parameter in section 2. Section 3 will present some example spectral energy 
distributions (SED) showing the significant differences between the linear and 
the nonlinear approach. In section 4 we will give examples of the inclusion of 
external photon fields into the cooling scenario, again altering the results. We 
will conclude in section 5. 



2 The injection parameter 

Schlickeiser, Bottcher & Menzler (2010) used the cooling term Q and a source 
term 5(7, t) — qo S (7 — 70) 5 (t) to solve the differential equation 0. As stated 
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above, the solution strongly depends on the injection parameter a, giving for 
a < 1 

n( 7 , t, a « 1) = g0 H [70 - 7] * (7 - 1 + ^ J , (3) 

which is solely determined by the linear synchrotron cooling. 

For a S> 1 the solution is divided into two parts for early and late times. 
For early times one finds 

n(7,7o J t<«c,a>l) = g -H'[7o-7] H[t c -t] 5 (7 



(1 + 3A go7o 3 1/3 



(4) 

reflecting the nonlinear cooling. For later times 

7S \ 



n(7,7o,* > tc,a> 1) = qo H[~/ B -7] H[t-t c ] 5 7- -j-^-g ,. 

(5) 

which is a modified linear solution. The critical time i c = (37b£>o) \ and 
7s = 7o/a- 

We call a the injection parameter, since it is defined as the square-root of 
the ratio of the non-linear cooling term to the linear cooling term at time of 
injection, i.e. 

2 \%sc{t = 0)\ A q j$ 

a = = ^0 (6) 

For a < 1 we see that the linear cooling dominates, resulting in the solution 
@. If a > 1, the non-linear cooling at least initially dominates, giving the 
solutions Q and ([5| for early and late times, respectively. The time t c marks 
the transition from non-linear to linear cooling. 

Since the Compton dominance, i.e. the ratio of the peak values of the two 
components in the SED, should be related to the ratio of the cooling terms, we 
expect a dominating synchrotron component for a 1, while for a 3> 1 the 
SSC component should show higher fluxes. 



3 Example SEDs 

In figure [l] we present two example SEDs f'(v'). They were calculated by 
Zacharias & Schlickeiser (2012a) by first calculating the intensity spectrum for 
each case and then integrating these spectra with respect to time, giving the 
fiuence spectra. The SED follows by multiplying the fluence with the frequency 
v. With the help of the Doppler factor S we transformed our spectra to the 
system of rest of the host galaxy (primed quantities) . 

The left panel in figure [T] shows the case for a -C 1, and one can clearly see 
that the synchrotron component dominates the SSC component. 

The right panel shows the case for a ^ 1, and here the SSC component 
dominates the synchrotron component as expected. Secondly, both components 
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Figure 1: Left: SED for a < 1 in the Thomson-limit. Right: SED for a > 1 
in the Thomson- limit. Parameters for each case are given at the top. Indicated 
are the Compton dominance (vertical line) and ratio of the peak frequencies 
(horizontal line), as well as the power-laws. 



exhibit broken power-laws, which are a natural result of the nonlinear cooling. 
We note the important fact that these breaks are due to the cooling term and do 
not depend on the underlying electron distribution (c.f. Zacharias & Schlickeiser 
2010 for the synchrotron component). Thus, we can explain the breaks in the 
SEDs without the need for fancy electron distributions, which are often given 
without any theoretical justification (e.g. Abdo et al. 2011). 

4 External photon fields 

Especially flat spectrum radio quasars (FSRQs) show strong thermal compo- 
nents, albeit the beamed nonthermal radiation. Thus, the electrons of these 
sources could also be cooled by inverse Compton cooling with these thermal 
components, like photons directly from the accretion disk (Dernier & Schlick- 
eiser 1993), from the broad line regions (Sikora et al. 1994) or the dusty torus 
(Blazejowski et al. 2000, Arbeiter et al. 2002). 

Zacharias & Schlickeiser (2012b) included this type of inverse Compton cool- 
ing into the above given model by adding the inverse Compton cooling term 

(Ted « 4.4 x 10- 8 ( <c \ Ilfs- 1 , (7) 
\ erg cm 6 J 

to the cooling rate Here u' ec is the energy density of the external photon 
field in the galactic frame, and Tf, is the Lorentz factor of the plasma blob in 
the jet. 

After calculating the resulting external Compton intensity and fluence, Zacharias 
& Schlickeiser (2012b) were able to present the complete SEDs, of which two 
examples are given in figure [2] 

The plots depend now also on the relative strength of the linear cooling terms 
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log(v') log(v') 



Figure 2: Left: SED for a < 1 and l ec < 1. SED for a > 1 and ^ ec > 1. 

Parameters for each case are given at the top. 

lec = Ijecl/ljsynl, and the normalized external photon energy e ec (assuming a 
monochromatic photon source). 

The main result is that the final SED depends strongly on the free parame- 
ters, and that both a dominating SSC or EC component are possible. Zacharias 
& Schlickeiser (2012b) argue, therefore, that the SSC component should not be 
neglected in blazar SEDs, and the correct time-dependent treatment should be 
applied. 

5 Conclusions 

In this proceeding we have given the results of recent analytical investigations 
concerning the time-dependent nature of the SSC process and the consequences 
for blazar SED modeling attempts. We first presented the results of Zacharias & 
Schlickeiser (2012a), in which we calculated the SED including the synchrotron 
and the SSC component. Then we proceeded with the recent work of Zacharias 

6 Schlickeiser (2012b), where we added the effect of external Compton cooling, 
yielding a third broad component in the SED. 

To conclude, we can say that our investigation has clearly shown the impor- 
tance of the inclusion of the time-dependent nature of the SSC cooling term, 
especially in blazar flaring states, which are far from steady-state. The resulting 
spectra differ significantly from the usual linear approach, where these effects 
are not taken into account. Secondly, due to the broad range of free parameters 
in the EC model, the SSC component should not be dismissed beforehand while 
modeling blazars. As we were able to show here, the SSC component might be 
as important as the EC component, if the correct time-dependent treatment is 
applied, which is especially necessary in flaring states. 

We aim to calculate the light-curves and the optical depths in a future work 
in order to further highlight the importance of the time-dependent SSC cooling 
treatment. 
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